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Abstract 

I present a simple analytical model describing the normal state of a super- 
conductor with a pseudogap in the density of states, such as in underdoped 
cuprates. In nearly two-dimensional systems, where the superconducting tran- 
sition temperature is reduced from the mean-field BCS value. Cooper pairs 
may be present as slow fluctuations of the BCS pairing field. Using the self- 
consistent T-matrix (fluctuation exchange) approach I find that the fermion 
spectral weight exhibits two BCS-like peaks, broadened by fluctuations of 
the pairing field amplitude. The density of states becomes suppressed near 
the Fermi energy, which allows for long-lived low-energy Cooper pairs that 
propagate as a sound-like mode with a mass. A self-consistency requirement, 
linking the width of the pseudogap to the intensity of the pairing field, deter- 
mines the pair condensation temperature. In nearly two-dimensional systems, 
it is proportional to the degeneracy temperature of the fermions, with a small 
prefactor that vanishes in two dimensions. 
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I. INTRODUCTION 



One of the strange features of high-temperature cuprate superconductors, a normal-state 
gap in the density of states, remains a subject of controversy. A great deal of experimen- 
tal evidence (NMR,§ optical conductivity,! specific heat, I and, most recently, angle- resolved 
photoemissionSli) seems to indicate that the superconducting energy gap survives the tran- 
sition to the normal state and disappears only at a considerably higher temperature. Even 
the anisotropic character of the gap is preserved. This behavior is characteristic of "under- 
doped" compounds, in which doping of additional hole carriers increases the temperature of 
the superconducting transition Tc. There is no generally accepted theoretical model of this 
phenomenon, although arguably the most popular tentative explanation puts the blame on 
superconducting fluctuations near the transition point, which are expected to be enhanced 
in these highly anisotropic, almost two-dimensional, materials. While Monte-Carlo simula- 
tions of the attractive Hubbard model in two dimensions provide numerical evidence! for the 
pairs-above-Tc scenario, no analytical description of the pseudogap regime has been offered 
so far. The goal of this paper is to provide such a sketch. 

Suppose that the pseudogap in cuprates does have the same origin as the superconducting 
gap — the scattering between fermion and hole states with the charge ±2 released in the form 
of a low-energy Cooper pair. This picture rests on certain assumptions. Simply saying that 
Cooper pairs are present above Tc as a propagating mode is not enough. It is important that 
the fluctuating pairing field A look frozen on the time scale 1/|A|, otherwise the pseudogap 
will not be formed. Since the typical frequency of free propagating bosons in thermal 
equilibrium is given by their temperature, the pseudogap is expected to exist at temperatures 
T ^ To, where Tq is loosely defined to be of order \A\/kB- In general, it may also be 
necessary to require the spatial coherence of the pairing field over a pair size ^ v/\A\. In 
this particular model, however, the boson velocity is of order of the Fermi velocity v, so that 
as long as the pairing fluctuations are slow in time, they are also slow in space. 

In a somewhat arbitrary way, Tq can be identified with the transition temperature calcu- 
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lated in the BCS theory. As in any other mean-field theory, the onset of a long-range order 
coincides with the appearance of a short-range order (formation of pairs), the assumption 
being that the phase of the order parameter will lock up as soon as there is a non-zero 
amplitude. This, of course, breaks down in two spatial dimensions, as shown long ago by 
Hohenberg.i In a highly anisotropic, almost two-dimensional system, long-range order sets 
in at a temperature Tc ^ Tq. Above Tq, Cooper pairs decay; between Tq and Tc, they 
represent a propagating mode; finally, below Tc, they form a Bose condensate. In the inter- 
mediate range, at least when Tc < T <^ Tq, they represent a slowly fluctuating BCS pairing 
field and thus create a pseudogap. A well-defined pseudogap regime may be a peculiarity of 
low-dimensional systems. 

A theoretical framework for treating the interplay between fermions and their boson-like 
bound states is thus needed. Following pioneering works by Eagles Jl Leggett,! and Nozieres 
and Schmitt-Rinkjli several approximate methods have been suggested: the boson-fermion 
model,00 functional integration,lii"Ei and the self-consistent T-matrix approximation.EUll 
In most cases, however, one has to resort to numerical computations. The present work is 
written with the purpose to provide an analytical sketch of a normal state with a pseudogap. 
As such, it inevitably contains further simplifications, which hopefully do not alter the nature 
of the problem: (1) The decay of low-energy bosons is neglected. (2) A clear separation of 
energy scales is assumed: 

fcsTc < |A| < ep (1) 

[ep is the Fermi energy). 

Using the self-consistent T-matrix approach, I derive approximate propagators for 
fermions and Cooper pairs. The width of the pseudogap is determined by the mean-square 
fiuctuation of the pairing field {\A{x,x')\'^) in a thermal ensemble. This expectation value 
depends, among other things, on the energy spectrum of Cooper pairs, which in turn is a 
function of the fermion energy spectrum. A closed set of coupled equations results, allowing 
one to determine the condensation temperature of pairs. The so determined Tc is propor- 



tional to the fermion density and inverse mass, as seen on the Uemura0 plot. This is despite 
the fact that, as long as |A| ^ ep, the system is not in the limit of local (tightly bound) 
pairs. I also give a reason why quasiparticle peaks are quite broad near the Fermi surface, 
even when thermal fluctuations are slow {ksTc ^ If correct, this sketch may provide 

a simple way to understand two puzzling features of underdoped cuprate superconductors - 
the pseudogap and the doping dependence of Tc - from a unified standpoint. 

The paper is organized as follows. The conserving T-matrix approximation is outlined in 
Section |T|; a bosonic propagator for a Cooper pair is introduced for the case of a "separable" 



interaction vertex. Section 111 contains a derivation of the fermion propagator in the presence 



of a slowly fiuctuating pairing field. An exact model of Section |^ reveals an important 
difference between a true superconducting gap and a pseudogap, stemming from the quantum 
nature of pairing fluctuations. Low-energy properties of Cooper pairs in the presence of a 
pseudogap and their condensation temperature are derived in Section 0. The effect of fast 
fluctuations is estimated in the Appendix. 

II. SELF-CONSISTENT T-MATRIX APPROXIMATION. 



A good description of the self-consistent T-matrix approximation can be found in Ref. 16 



The relation to other approximate methods is outlined in Ref. |T8. 

The conserving T-matrix approximation is somewhat similar to the well-known Hartree- 
Fock principleJni The latter neglects any correlations between interacting particles (except 
for statistical ones) and describes the motion of independent entities in a self-consistent 
potential. The T-matrix approach includes pairwise correlations between colliding particles, 
which are particularly important when two-particle bound states are formed. The form of 
these correlations is inferred from the exact solution of a similar problem with two particles 
in vacuum, when the two-particle Green's function 

G^ix, x'- y\ y) ^ -i{T[i^{x)i^{x')i^\y')i^\y)]) (2) 

can be expressed in terms of the one-particle Green's function G{x;y) = —i{T[ip{x)i/j'^ (y)]) 



and the two-body T matrix [Fig. |T](a)]. In the many-body case, the latter is defined as a 
solution of the Bethe-Salpeter equation [Fig. |l|(b)]: 

T{P\k, k') = V{P\k, k') + ^ Y.T{P% k")G{P/2 + k")G{P/2 - k")V{P\k\ k'). (3) 

Here P = {Q, P) is the total 4-momentum of two fermions, k and k' are relative 4-momenta, 
and V is the four-dimensional volume (with an imaginary time dimension). By reducing G2 
to a functional of G, one breaks the infinite hierarchy of equations for n-particle Green's 
functions because the resulting approximate Dyson equation G = G^^^ + GT,G^^^ contains 
no higher-order Green's functions [Fig. |](c)]: 

m = E nP\P/2 - k, P/2 - k)G{P - k), (4) 

K p 

Equations § and ^ form a closed set of equations that could be solved, at least in principle.lli 
The use of full fermion propagators (instead of bare ones) in these equations is the 
main difference from the earlier approach of Thouless,0 equivalent to the BCS theory. In 
addition to having conservation laws built in,ii the self-consistent approach is a step towards 
including the influence of Gaussian fluctuations on the phase transition. Near the transition 
temperature, fermions start to feel the presence of the emerging bosonic excitations, which 
in turn influences the process of pair formation. Dressing of the fermionic propagators in 
Eq. ^ allows one to account for this feedback, at least to some extent. The use of the 
dressed fermion propagator in the self-energy equation (^ increases the number of terms in 
the perturbation series for G{k), which turns out to be important for treating properly the 
quantum nature of fluctuating Cooper pairs. 

A. Propagator of a Cooper pair. 

Although the T matrix depends on three independent momenta, the non-trivial depen- 
dence should be associated with the conserved total momentum P. It is convenient to write 
the interaction vertex formally as [Fig. |I|(d)] 
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V{P\k,k') = u\k)D^^\P)u{k') = / u*{k\x)D^^\P\x,x')u{k'\x') d^xd^x', 



(5) 



where x and x' are some relative coordinates. One well-known example of such factorization 
is the case of a separable instanteneous potential V{P, k, k') = gv(k)v(k'). This formalism 
is particularly convenient for lattice models with finite-range instanteneous interactions. In 



the case of the t-J model, for instance, 



V{P,k,k') = (^1 cosA;^ cosfcy) 



+00 
-J 
-J 



1 

cos k' 



cos k' 



(6) 



in the singlet channel and zero in the triplet one. The infinite on-site repulsion imposes the 
constraint of no double occupancy. Phonon-mediated attraction can also be written in the 
form (H) with u{k\x) = e*'^'''""^*. 

The T matrix can now be found using the ansatzT{P\k, k') = u\k)D{P)u{k') [Fig. 0(e)]. 
The new matrix D{P\x^x') satisfies the Dyson equation D = D^^^ + DUD^^^ with the 
polarization matrix 



n(P) = y E u{kW{k)G{P/2 + k)G{P/2 - k). 



(7) 



Clearly, D{P\x,x'), which contains all the non-trivial information about the T matrix, can 
be regarded as the bosonic propagator of a Cooper pair whose internal structure is described 
by the dependence on the relative coordinates x and x' . 

The factorization procedure described above allows one to write the Bethe-Salpeter equa- 
tion for the T matrix in the form reminiscent of the Dyson equation for a doubly charged 
bosonic particle — see Fig. |l](f) — without introducing spurious degrees of freedom. The 
Thouless criterion for pair condensation, T(P) = cxo at P = 0, simply means that the 
quasiparticles with P = have zero excitation energy. 
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III. FERMION PROPAGATOR IN THE T-MATRIX APPROXIMATION. 



The fermion self-energy @ can now be rewritten using the boson propagator D{P) 
[Fig. ig)]: 

S(fc) = --Y. u\P/2 - k)D{P)u{P/2 - k)G{P - k). (8) 
V p 

This diagram illustrates two important points about the nature of the T-matrix approxi- 
mation. (1) The irreversible decay of a fermionic quasiparticle is determined by the density 
of states "hole + Cooper pair" at a given 4-momentum. Note that both outgoing lines are 
dressed and thus represent actual excitations of the system. This should be contrasted to 
the functional-integration approach at the usual one-loop leveliil where the fermion self- 
energy is expressed in terms of the hare fermion propagator. We will return to this point 
later in Sec. |rV|. (2) For a slowly fluctuating pairing field, the sum in Eq. ^ is dominated by 
the region near P = 0, so that a fermion with 4-momentum k is coupled mostly to the hole 
with the same 4-momentum. By replacing G{P — k) with G{—k), we obtain 

m « -|A(.)|=G(-.) . (9) 

where 

\m\' ^Y;^ n\Pl2 - k)D{P)u{P/2 - k) (10) 
p 

The approximation made here may appear rather crude. Essentially, the scattering of a 
fermion into a continuum of hole states is replaced by the coupling to a single hole state 
with the same 4-momentum, a situation reminiscent of the BCS superconductor. If a bare 
propagator were used for the hole, the resulting Bogoliubov quasiparticles would be stable. 
However, the use of a full propagator G{~k) already makes the lifetime finite. In this 
case, neglecting temporal and spatial fluctuations of the pairing field does not look so bad, 
especially in a low-dimensional system. 

The approximate equation for the self-energy can be readily solved.^ Iterating it once 
results in a quadratic equation for S(A;). Since |A(A;)p = |A(-A;)|2, the dressed propagator 
is 



G(k) = 

I + J1 + 4|A(A;)|2G(o)(A;)G'(o)(-A;) 



The branch of the square root is fixed by the requirement that G{k) G^'^\k) as |A(A;)| 0. 

In the simplest case of a short-range instanteneous attractive potential, A(A;) can be 
replaced by a constant and we have 



1 I 

G{u,k) 



1 lAP 



-1 



uj^ — et 



(12) 



where ek is the bare fermion energy. It is readily seen that the spectral weight of the dressed 
fermion is distributed over two finite regions < u;^ < + 4| Ap as follows: 

+ Ckl 



^(u;,k) 



27r|A|2 \ 

This distribution is reminiscent of the smeared BCS peaks at u;^ = + |Ap. In fact, the 



ratio of the spectral weights (p^) at = iyek + I^P the same as in the BCS case. 

Since the spectral weight is pushed away from the Fermi surface, a pseudogap opens up 
in the total density of states Af{uj) (Fig. 0): 

Afiuj)=Mofiiu/2\A\), (14) 

where A/q is the density of states in the free case and the function /(x) can be expressed in 
terms of complete elliptic integrals: 

{(4/7r)a;E(a;) for < x < 1, 
(15) 
(4/7r)[K(l/x) - D(l/x)] for x > 1. 

It vanishes linearly for small values of x and approaches 1 as x 00. 

IV. IDEAL FERMIONS INTERACTING WITH AN INDEPENDENT PAIRING 

FIELD. 



The T-matrix approximation amounts to summing up an infinite number of terms of 
the perturbation series. Yet an even "greater" infinity of terms is left out. In such a case, 
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when approximations are not easily controlled, it is desirable to check the result using some 
exactly solvable model. 

Consider a free fermion gas interacting with a classical external pairing field, which is 
represented by the action term 



The field A(x, x') can be interpreted as the wavefunction of a Cooper pair with its constituent 
fermions at x and x'; the interaction term above removes two free fermions and forms a bound 
state. 

The U{1) symmetry, related to the conservation of charge, makes the phase of the com- 
plex field A unobservable. When pair formation lowers the energy of the system, the density 
of pairs |Ap may become large. There are two distinct possibilities in this case. The first 
one occurs when the U{1) symmetry is sponateously broken and A chooses some direction 
in the complex plane with the well-known results (the BCS superconductor with a gap and 
a Bogoliubov sound mode). The second possibility is a symmetric phase with a fiuctuating 
field A. The quantum nature of these fiuctuations leads to a non-trivial effect: a significant 
broadening of the fermion spectral weight. 

Quantum fluctuations in this model can be implemented by using an ensemble of pairing 
fields A with a symmetric distribution. Although every single measurement of A produces 
a deflnite non-zero result, (A) = upon averaging over a long series of such measurements. 
This procedure is Feynman's way of quantizing a classical system.!! A closely related exactly 
solvable model for the Peierls gap in quasi-one-dimensional metals has been discussed by 



A generic diagram of the perturbation series for the fermion propagator contains n incom- 
ing and n outgoing dashed lines representing the pairing fleld [Fig. ^(a)]. If {A{x,x')) ^ 0, 




(16) 




following the analytic solution of Sadovskii.l 



A. Broken phase. 
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there is also the anomalous fermion propagator, which gives rise to the Meissner effect 
and superconductivity. Consider the extreme case when there are no fluctuations in the 
ensemble: A{x,x') = {A{x,x')). The dressed fermion propagator is simply a geometri- 
cal series with a single term in each order. For an instanteneous short-range pairing field 
{A[x,x')) = A6{x — x') (uniform for simplicity), 

2 



G{k) 



E 



A 



(17) 



Although the geometrical series J2^=oZ^ diverges when \z\ > 1, its analytic continuation, 
the function (1 — z)~^, has only one singular point z = 1. The BCS result is recovered: 



G{k) 



A 



The divergence of the perturbation series (|T^ in the range < u;^ < + 



A 



(18) 



is associated 



with the transfer of the pole of the propagator from = ek to the new locations at cu^ 

2 



4 + 



A 



B. Symmetric phase with quantum fluctuations. 

The structure of the perturbation theory changes considerably in the symmetric phase. 
Assuming that the statistical distribution of the pairing fields is Gaussian, we can write any 
average that involves a product of n fields A{xm,x'^ = Am and n fields A*{ym,y'm) = 
as a sum over all possible products of pairwise averages: 

(AtAiA^As . . . A:A„) = 5:(AtApi)(A;Ap2) . . . (A:Ap„), (19) 

p 

where P represents a permutation of the numbers 1 . . .n. Upon introducing a dashed line 
for each pair average (A*Aj), we find n\ different diagrams with n such lines [Fig. ^(b) and 
(c)]. Again, take the example of a uniform instantaneous short-range pairing field: 

{A*iy,y')Aix,x')) =WSix - x')6iy ~ y'). (20) 
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Each single field breaks the phase symmetry and produces the effect described above. Av- 
eraging over the ensemble, however, restores this symmetry; the perturbation series now 
contains n! identical terms in the n-th order: 



1 ~ / IA|2 \ 

^(^) = — E^M-^ • (21) 



Clearly, the result is different from (^7^. Expression (|2lD was obtained in Ref. ^ in connec- 
tion with the pseudogap in the boson-fermion model - as a leading term in two dimensions. 
Although the series 



oo 



J2 nlz"" (22) 

n=0 

diverges for any z 7^ 0, it represents an asymptotic Taylor expansion of the function 

The divergence of the series is related to the nonanalyticity of f{z) on the positive side 
of the real axis, where it has a cut. Rather then working with a series that converges 
only in the trivial case z = 0, we can sum up the diagrams for a single uniform field A 



- as in Sec. [IV A| - and then average the result over a Gaussian ensemble with the weight 



exp -|A|V|A|2WA*dA: 



G{k) = — I — , (24) 



k 



which indeed has (]2l|) as an asymptotic expansion. The quasiparticle spectral weight, con- 



centrated at = ek in the free case, is now spilled into the region uj"^ > e\ (Fig. 

^(-,k) = ^^exp(-^^), c.^>4. (25) 

The density of states exhibits a pseudogap (Fig. somewhat more pronounced than it 
was found in the T-matrix approximation: it now vanishes as u"^ near the Fermi level. 

The results of this section can be understood in a simple intuitive way. Even though 
we have suppressed spatial and temporal fluctuations of the pairing field, its strength does 
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not yet have a definite value; the fermion pole, instead of being shifted to a definite new 
position, is smeared into a cut along some part of the real axis. Nevertheless, the existence 
of a scale |Ap in the field distribution results in the fermion spectral function reminiscent 
of the BCS one. 



C. Relation to the T- matrix approximation. 

The T-matrix equation for the fermion self-energy (|^), which includes a dress ec? fermion 
propagator, is equivalent to partial summation of the diagrams in Sec. |1V B\ only those 
without intersecting dashed lines are included. For instance, diagram (b) in Fig. |^ is present, 
while (c) is left out. The number of diagrams in the n-th order satisfies the recursion relation 

n 

Cn+1 = ^ CmCn-m (26) 
m=0 

(Co = 1 by definition) with the solution 

C'n=n 4-— -, (27) 



m=l 



171+ 1 



which grows as 4" for a large n. Therefore, the series 

oo 

E Cn^" (28) 

n=0 

has the radius of convergence equal to 1/4 and the Green's function 



1 °° / IA|2 \ 
G(A:) = -— (29) 



is analytical (hence real) when u"^ > + 4|Ap. The spectral weight is contained in the 



finite region < < + 4| Ap. 

The above discussion shows that the T-matrix approximation includes, at least partially, 
effects of quantum fluctuations of the pairing field on the fermion propagator. This is 
achieved through the use of a dressed fermion propagator in the equation for the fermion 
self-energy. The number of terms in the resulting perturbation series grows quickly enough 
to make the result non-analytical in a larger region than in the BCS broken phase. Inclusion 
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of a greater number of terms (the full series) only slightly enhances this effect. This should be 
contrasted to the situation when a bare fermion propagator is used in the fermion self-energy 
(y). Since there is only one graph in each order of the perturbation series for the propagator 
G{k), such a treatment essentially amounts to throwing away the quantum nature of Cooper 
pairs. No broadening is then found in the limit of slow fluctuations. 



V. COOPER PAIRS IN A PSEUDOGAP. 

A. The propagator of a Cooper pair. 

With much of the spectral weight removed from the vicinity of the Fermi surface, low- 
energy bosonic excitations may represent a propagating mode. In this respect, such a normal 
state resembles a BCS superconductor just below Tc, where the opening of a true gap inhibits 
the decay of pairs.0 Therefore, a reasonable estimate of the boson energy spectrum can be 
obtained from Eq. ^ by using the dressed fermion propagator in the BCS form (|18|), which 
constitutes the two-pole ansatz of Ref. 0. Inclusion of fermion lifetime effects makes Cooper 
pairs unstable but is not expected to produce qualitative changes in the energy spectrum. For 
a short-range instanteneous attraction, D^^\P) = g < 0, 

D-\P) = + 1 X: G{k)G{P - k), (30) 

with 

G{k)^G{u;,k)= ""t'W - (31) 
Lu^ — e^ — \A\^ 

At this level of approximation (no boson decay because of a full fermion gap), D^^{Q, P), 
analytically continued from Matsubara frequencies to the rest of the Q plane, can be ex- 
panded in powers of Q and P: 

fl-,aP) ^ fl-'(0, 0) . 'Jg^n^ + «)P= + . . . (32) 

or 
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The mass term appears because the Thouless criterion no longer apphes: D^^{0,0) 7^ in 
the normal state. The existence of two Cooper-pair poles at a given momentum P, at the 
frequencies Q = ±V M^s"^ + P^s^ = ±i?p, reflects the fact that correlated propagation is 
possible not only for two fermionic particles, but also for two holes. If a pair is made out of 
fermionic excitations just around the Fermi surface, where the density of states is the same 
for particles and holes, the residues of the two poles are equal (up to a sign). A varying 
density of states will make the two poles less symmetric, both in terms of their residues and 
positions, but still low-energy Cooper pairs will have two poles. 

This trend can be clearly seen in the numerically obtained two-particle density of states 
plotted in Fig. 4 of Ref. Low-momentum pairs have two somewhat asymmetric peaks 
dispersing with momentum. The asymmetry is related to the instanteneous nature of the 
interaction in the attractive Hubbard model: the number of particle states involved in the 
formation of a pair greatly exceeds that of hole states at low fermion densities. No holes 
contribute to the formation of pairs with higher momenta, so that there is no pole at negative 
frequencies in this case. Such pairs, however, represent fast fluctuations of the pairing field 
and are irrelevant to the formation of the pseudogap. The numerical data also show an 
increasing lifetime of Cooper pairs as P ^ 0, reflecting a stronger depletion of the fermion 
density of states near the Fermi energy. 

When fermionic excitations across the gap are frozen out, the Cooper-pair propagator 
at zero momentum is 

B-.(n.o)^9-'-/>W*94^ (34) 

where e = \J + With the exception of the zeroth-order term D^^(0, 0) = —Z~^M'^s'^, 
coefficients of the Taylor series (in powers of Vt) are insensitive to large-e behavior of the 
density of states, and are thus determined by the lower energy scale, i.e., the gap width |A|, 
provided that variations of the bare density of states are small on that energy scale. In this 
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weak-coupling limit, the terms odd in Q vanish because M{e)e is an odd function near e = 0. 
Quartic and higher-order terms can be dropped as long as Vl? <^ |Ap. Then 

D-\Q,^)^ Z-\n'' - M\s^), (35) 



where Z ^ is of order Ao/IAp. More exactly, 



Z = 12|A|VAro. (36) 

Expansion in powers of P allows one to determine the characteristic speed of bosons 
s, which is expected to be smaller than the Fermi velocity v on physical grounds. In the 
weak-coupling limit, we recover the result of Bogoliubov, = v'^/d in d dimensions. This is 
roughly consistent with the dispersion of Cooper-pair poles at low momenta inferred from 



the data of Ref. IS . 



At a given temperature, two as yet undetermined low-energy parameters of the boson 
field - the mass M and the pole residue Z - control the intensity of fermion scattering by 
the pairing field and thus determine the (also unknown) fermion energy gap, as discussed in 
some detail below. In addition, Eq. relates Z to the energy gap. Using this information, it 
should be possible to determine two out of the three unknown parameters. In the following 
section, the mass M is found as a function of temperature. 



B. Evaluation of the fermion self-energy 

We now return to the analysis of the fermion self-energy (H) using the approximate boson 
propagator (|33|). The following expression for the fermion self-energy results: 

The sum over the bosonic Matsubara frequencies Q has been converted into an integral 
around the poles of (e^*^ — 1)^^; Ep = ^ M'^s^ + P^s^ > is the boson energy. It is 
convenient to write the fermion propagator in terms of its spectral weight: 
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Gi^,k) ^ r (38) 

J —oo UJ — 6 

After deforming the integration contour in the standard way, one finds the self-energy as a 
sum of two terms, which correspond to emission (+) and absorption (— ) of a Cooper pair: 

The part proportional to the boson occupation number N{E-p) can be regarded as induced 
emission or absorption with the rate proportional to the spectral intensity of the pairing 
field N{Ep)/Ep. The fermionic term n(±e') is responsible for "spontaneous" processes. 

If the pseudogap formation is governed by low-energy modes, the induced term is domi- 
nant since N{Ep) ^ n{±e') at low energies. The rate of spontaneous emission and absorp- 
tion of pairs is estimated in the Appendix; it can be neglected in c? = 2 dimensions. Upon 
shifting the integration variable e' =F Ep e', 

Bosons are restricted to low-energy modes with Ep of order or less than ksT . In the 
case when the fermion spectral weight is distributed over a wider interval of energies, the 
variation of the fermion spectral weight on the energy scale of fc^T can be neglected. For 
the same reason, we can neglect the variation of the fermion momentum in the integrand, 
provided that fermion and boson velocities are of the same order (recall that = v'^/d). 
The integration over e' yields the fermion propagator G{—uj, — k) and the self-energy has the 
form inferred previously (|^): 

sL1(u;,k) + sL-i(^,k) ^ -WG{-uj, -k), (41) 

where the average fluctuation (intensity) of the pairing field is 



|A|2 = Z / , ^ (42) 

' ' J {2nY Ep ^ ' 

Essentially the same result can be obtained by considering a classical field A(t, r) with 
the Lagrangian density determined by (^31), 
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£(t,r) = Z-i |5A(t,r)/5tp- |VA(t,r)|2-Ms2|A(t,r)|2 , 



(43) 



in thermal equilibrium. Application of the equipartition theorem yields 




(44) 



which is the classical analogue of Eq. ^ 



C. Condensation temperature. 



By substituting the value of the boson residue (|36D into Eq. we obtain the self- 
consistency condition mentioned above. Remarkably, the width of the pseudogap cancels out 
and the resulting equation implicitly determines the boson mass as a function of temperature: 



where Ep = \J M'^s^ + P^s^, N(E) = (e^^ — 1)^^ is the boson occupation number, and 
A/q is the density of bare fermion states. It is immediately obvious that there is no Bose 
condensation in ci = 2 dimensions: for M = 0, the integral diverges at the lower limit. 
In d = 3 dimensions, this equation predicts a condensation temperature fcsTc of order ep- 
This result should not be taken at face value because it was derived for a weakly attractive 
degenerate fermion gas. (Also, as discussed in the Appendix, the absence of a low energy 
scale implies that pairing fluctuations are not slow.) Nevertheless, the existence of such a 
high temperature scale is justified if one considers the limits of weak and strong attraction 
between fermions. In the weak-coupling limit, there are no pairs at the Fermi temperature, 
they form at a much lower temperature fc^To ep; therefore, the long-overdue condensation 
occurs immediately, which explains why the mean-field BCS approach works so well. In the 
opposite limit of local pairs (turned into hard-core bosons) , condensation indeed occurs when 
the gas becomes degenerate. 

Next, we estimate the condensation temperature in d = 2 + e dimensions, which may be 
relevant to highly anisotropic cuprate superconductors. Recalling the asymptotic behavior 
of the Riemann zeta-function ^(1 + e) ~ 1/e we obtain 



/ 



{27iY Ef ~ T2 ' 



(45) 
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kBT, ~ eep/Q (46) 

when £ ^ 1. For comparison, the condensation temperature of a Bose gas made of very 
small fermionic pairs is only higher by a factor 3/2. (In both cases, ep is the Fermi energy 
of the ideal fermion gas with the same mass and density of particles.!!) It is thus clear 
that the condensation of Cooper pairs in the presence of a pseudogap occurs well below the 
fermion degeneracy temperature and a clear separation of energy scales ([I|) is possible for 
moderately weak attraction strengths. If the attraction is strong, we end up with tightly 
bound pairs, in the local boson limit. In the case of a very weak attraction, the BCS pair 
formation temperature Tq is lower than Tc of Eq. ^ the BCS model takes over. 

VI. CONCLUSION 

I have considered the interaction between fermions and Gaussian pairing fluctuations 
(Cooper pairs without self- interaction). A well-defined pseudogap regime is found ind = 2+e 
spatial dimensions, when the condensation temperature is much lower than the BCS mean- 
field one. The quantum character of Cooper pairs, related to the unobservable nature of the 
pair wavefunction A(x, x'), leads to considerable broadening of the Bogoliubov quasiparticle 
peaks near the Fermi surface. Even when non-uniform configurations of the field A{x,x') 
are suppressed at low temperatures, fluctuations of its amplitude blur the strength of the 
pairing field, making the position of the quasiparticle poles at 07^ = 6^ + |A(/c)p uncertain. 

A well-pronounced suppression of the fermion density of states near zero energy makes 
room for long-lived pair states inside the gap. At low momenta and frequencies, they rep- 
resent a Bogoliubov sound-like mode with a non-zero mass. The mass is determined from a 
self-consistency condition that links the width of the pseudogap to the mean fluctuation of 
the pairing amplitude. Zero mass signals the onset of Bose condensation. Ind = 2+e dimen- 
sions, the condensation temperature scales as eep/Qi which is 2/3 times the condensation 
temperature for the corresponding ideal Bose gas. This may explain the doping dependence 
of Tc in underdoped cuprates. 
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Approximations made in this work are twofold. First, an infinite hierarchy of equations 
for n-particle Green's functions is replaced with the self-consistent T-matrix approach of 
Brueckner, which partially accounts for Gaussian pairing fluctuations. By solving an exact 
model with an independent pairing field, I show that the self-consistent T-matrix approxi- 
mation takes into account the amplitude fluctuations of Cooper pairs in the normal state, 
a feature missing in other approximate models.lii0@ Further approximations made in this 
paper are necessary to obtain an analytic solution. Firstly, thermally excited non-uniform 
configurations of the pairing field have little effect on the fermion propagator, provided 
that a clear separation of energy scales, fc^Tc ^ |A|, exists. This appears to be the case in 
d = 2 + e dimensions, i.e., for highly anisotropic, almost two-dimensional systems. Secondly, 
a residual density of fermionic states in the pseudogap, leading to the decay of Cooper pairs, 
is assumed to have negligible influence on the pair energy spectrum. A numerical work 
using the self-consistent T-matrix approachlll seems to indicate that this is a reasonable 
assumption: low-energy Cooper pairs represent a propagating mode with two Cooper-pair 
poles dispersing with the pair momentum. In order to ascertain that propagating Cooper 
pairs in the pseudogap are not an artifact of the self-consistent T-matrix approximation, 
it is necessary to look for them in quantum Monte- Carlo simulations on two-dimensional 
lattices. 

There are many questions that remain open. For instance, in the case of a superconduct- 
ing order parameter with nodes, the density of fermion states is less strongly suppressed. 
Will Cooper pairs be stable enough to represent a propagating mode? Also, what causes 
the strong broadening of fermion quasiparticle peaks above the temperature at which pairs 
are formed? 
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APPENDIX: 

Spontaneous emission or absorption of Cooper pairs is not restricted to low boson mo- 
menta and thus represents the effect of fast fluctuations. The emission or absorption rate 
depends on the phase space available to the products of a fermion decay (a pair and a hole), 
which tends to be smaller in a lower number of dimensions d. This indeed is confirmed by 
a simple calculation below. In d — 3 dimensions, the imaginary part of Ssp(a;,k) is linear 
in cu and thus leads to a renormalization of the quasiparticle spectral weight. In d = 2 di- 
mensions, a constant imaginary part is introduced, which merely broadens the quasiparticle 
peak. Besides, the existence of a small parameter in the pseudogap regime (the ratio \A\/eF) 
makes this additional broadening insignificant in d = 2 + e dimensions. 

Let us estimate the impact of "vacuum" fluctuations on the self-energy of a fermion with 
momentum k and bare energy ek = 0. Inasmuch a hole created in the decay can be far away 
from the Fermi surface, we can neglect the influence of pair scattering on its spectrum and 
treate the hole as a stable particle with momentum — k' = k — P and energy 



k P IP 



2 



- = -ekH ^ « iPlvcos^ = -e'(|P|,^), (Al) 

m 2m 

neglecting the curvature of the energy surface ek; 6* is the angle between k and P. The boson 
energy is taken to be Ep — s\P\ = £'(|P|). Then 

Im S(^) u;, k = -Zn / - — A2 

^ ' ^ J {27cY 2E ^ ' 

Energy conservation requires that 

E=—^^. (A3) 

S-\-V COS 

For definiteness, let us assume that the frequency of the incoming fermion > 0. Then a 
pair can be emitted {E > 0) when cos^ > —s/v and absorbed [E < 0) when cos^ < —s/v. 
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Pair emission is accompanied by absorption of a second fermion with the energy in the range 
e' > —vuj/{v + s) = uj^. The second fermion can, thus, be taken from below the Fermi surface 
and the process will take place even at low temperatures! On the contrary, absorption of a 
pair requires emission of a second fermion with a negative energy e' < ~vuj/{v — s) = uj^. 
At low temperatures, such states are occupied and pair absorption is suppressed. 

In (i = 3 dimensions, the element of phase space reduces to E dE (ie'/(27r)^f s^. The 
integration over the boson energy E is elementary because a ^-function is present: 

Z 

The first term in brackets (emission of a pair) grows linearly with ut. By inserting Z from 
fBSl), we obtain 



Im Esp(u;,k) 



n{e')de'+ / n{-e')de' 



(A4) 



TT I ^ I '^UJ 

""^•-'"•''>~' 4(i + i/V3) 

as uj ^ +00. Thus, spontaneous pair emission leads to a significant renormalization of the 



fermion propagator (the ratio |Ap/e|, is not small for the pseudogap regime in = 3) and 
cannot be neglected. 

In d = 2 dimensions, pair emission has a smaller available phase space: 

Im S(+)(u;,k) = / ^ ' . (A6) 

The self-energy does not vary with uj if n(e') is replaced with the step-function. This channel 
of decay is characterized by the rate 



ImS3p(cu,k) = ^^, (A7) 



where sinha = ^1/2—1/2. Ind = 2 + e dimensions, |Ap ^ ep and this additional 
linewidth is small in comparison with the broadening caused by induced fluctuations, of 



order y |Ap. Therefore, the effect of fast spontaneous fluctuations can be neglected in 
d = 2 + e dimensions. 
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FIGURES 



FIG. 1. Diagrams of the self-consistent T-matrix approximation for the singlet channel. 

FIG. 2. Pseudogap in the fermion density of states in the limit of slow pairing fluctuations. 
Dashed line: T-matrix approximation. Solid line: exact model with Gaussian fluctuations. 

FIG. 3. Feynman diagrams for the exact model with a classical pairing field. Solid lines are 
bare fermion propagators, (a) A diagram of order n = 3. Dashed lines with circles are pairing 
fields A{x,x') and A*{y,y'). (b)-(c) Examples of diagrams generated from (a) after averaging over 
a Gaussian ensemble of pairing fields. Dashed lines are now pairwise averages {A{x,x')A*{y,y')). 

FIG. 4. Distribution of the fermion spectral weight A{u, k) in the model with Gaussian pairing 
fluctuations, ek is the bare energy of a fermion, |A| is the r.m.s. fluctuation of the pairing field. 

Also shown are positions of the fermion poles in a non-interacting fermion system (dotted line) 
and a BCS superconductor with the gap |A| (dashed line). 
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FIG. 1. 



25 



Mo 









1 1 1 1 1 1 1 






\\ / / \} 
^ \ \ / / \_ 




- 


\ \ / / 

\ \ / / 
M / ' 
M / / 
M /' 
M // 

A // 
A /' 

\ { 


- 




\ '/ 
V // 
\\ '/ 

\^ '/ 
\^ '/ 
\^ '/ 
\ ^ ' / 
\ \ / / 

\ ^ ' / 
\ \ / / 

1 1 1 VI'/ 1 1 1 




-2|A| 2|A| 



FIG. 2. 



26 



(a) ! ^ ! ^ 



(b) 



(c) >/'^ 
/ ■' / ^ 

— ^ — 6 — K^> 



< 



FIG. 3. 



27 




FIG. 4. 



28 



